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I . INTRODUCTION 


In  the  process  of  conducting  live  firings  of  missiles  equipped 
with  warheads  containing  numerous  bomblets,  a problem  arose  per- 
taining ^to  the  evaluation  of  the  effectiveness  of  the  bomblets  against 
various  targets.  The  targets  included  personnel  simulators  and 
surface-to-air  missiles. 

Briefly,  it  was  necessary  to  determine  the  probability  of 
kill  (PjJ  of  a group  of  target  items  when  some  or  all  of  them  were 

damaged  to  unequal  levels  of  kill.  In  essence,  it  was  required  to 
compute  the  probability  of  achieving  exactly  and/or  at  least  X 
kills  (successes)  out  of  a sample  of  N trials  when  the  individual  P,, 
values  were  not  necessarily  equal. 


II.  PROCEDURE 


Standard  textbook  methods1  use  the  following  equation  for 
computing  the  probability  of  X successes  out  of  N trials  when  the 
probability  of  success  p,  is  constant  from  trial  to  trial  and  the  N 
trials  are  independent: 


P(X;N,p) 


(1) 


whe re  X = 0,1,2  . . .N, 

M _ N! 

\\)  " X!  (N-X) ! 

and  q = 1 - p. 


However,  in  the  case  under  consideration  here,  the  values  for 
are  not  necessarily  equal  and  Eq.  (1)  cannot  be  used  to  compute  the 
probability  of  exactly  X successes  out  of  N trials,  P(X;N,p^)  which 


for  the  sake  of  simplicity,  we  will  express  as  P 


throughout  the 


P 


1 Freund,  J.E.,  "Mathematical  Statistics",  Englewood  Cliffs,  NJ : 
Prentice -Hall,  1962 
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remainder  of  this  report.  Instead,  consider  the  following  expression 
which  can  be  used  to  determine  ^or  unetIual  (°r  e(lual  > for  the 

trivial  case)  probability  values,  p^ : 


2N 

N 

-0-E 

J-l 

6i 

TTV1J 

i = l 

(2) 

where  the  6.  .'s  are  chosen  to 
i » J 

Eq.  (3)  for  each  j, 

equal  0 or  1 to  satisfy 

N 

m • £ 

<5. 

i.J 

2i_1 

(3) 

r“4 

II 

•H 

and  3.  is  defined  as 

J 

N 

3.  = 1 iff 

J 

£ 

i=l 

6.  . = X 
i.J 

(4) 

= 0,  otherwise. 


As  a simple  example,  consider  the  case  of  determining  the 
probability  of  at  least  1 success  out  of  3 events  and  exactly  1, 
exactly  2,  and  exactly  3 successes  out  of  3 events.  Using  Eq . (3) 

and  Eq.  (4)  where  j = 1 to  23  and  N = 3,  the  6.  ,'s  and  3.'s  are 
M i,3  J 

determined  as  follows: 

For  j * 1 , 

3 


j-‘  ■ £ 6i,j  2i_1 


i=l 


0 ' *1.1  2°  * S2,!  2‘  * *J,1  22 


In  order  for  the  above  equation  to  be  valid,  5 must  equal  0, 

A 9 * 


6 


6 must  equal  0,  and  <5  must  equal  0,  and  therefore  the 

2,1  *5 » 1 

3 

^ ^ 6 = 0.  Consequently  using  Eq.  (4),  for  X = 1 , 3j 


must 


i=l 

equal  0,  for  X = 2,  3^  must  equal  C and  for  X = 3,  must  equal  0. 
For  j = 2, 


j-1 

1 


3 


i*l 


♦ 


+ 


5 


2 


Similarly,  for  the  above  equation  to  be  valid,  6 must  equal  1 

I , z 


2,2 


must  equal  0 and 


6 . must  equal  0, 

O f Z 


and  therefore 


E 

i = l 


5. 


Again,  using  Eq.  (4)  for  X = 1 , 8?  must  equal  1,  for  X = 2,  8.,  must 
equal  0 and  for  X = 3,  3?  must  equal  0. 


For  j = 8, 

■i_i  ■ E *u  2‘'‘ 

i = 1 

7 - «1>8  2“  * *2,8  2‘  * *3,8 


Finally,  for  the  above  equation  to  be  valid, 


6,  0 must  equal  1, 
1,0 
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must  equal  1 and  6 


must  equal  1,  and  therefore  the 


2,8 


3,8 


£ 

6.  . = 3. 

Using  Eq. 

(4) 

» 

for  X = 1,  0 

must  equal  0, 

O 

for  : 

i=l 

0g  must 

equal  0 

and 

for 

X = 

3, 

e8 

must  equal 

1.  The 

formation  of 

6.  . 's 

from  Eq. 

(3) 

and 

the 

V 

s 

from  Eq.  (4) 

for  the 

preceding 

example 

in  which 

N = 

= 3, 

can 

be 

summarized  as 

follows : 

3 

j 

6l,j 

62,  j 

6 

3, 

V 6 

i i. 

6. 

1 1 

8. 

J 

B. 

1 

i=l 

X=1 

X=2 

X=3 

1 

0 

0 

0 

0 

0 

0 

0 

2 

1 

0 

0 

1 

1 

0 

0 

3 

0 

1 

0 

1 

1 

0 

0 

4 

1 

1 

0 

2 

0 

1 

0 

5 

0 

0 

1 

1 

1 

0 

0 

6 

1 

0 

1 

2 

0 

1 

0 

7 

0 

1 

1 

2 

0 

1 

0 

8 

1 

1 

1 

3 

0 

0 

1 

These  tabular  values  may  now  be  applied  to  Eq.  (2)  in  which  case 
the  probability  of  exactly  1 success  out  of  3 events  is  expressed  as: 


= 3 *(p°p°p°q1q1q1) 

\1/  1 123123 


e9*(p1p0p0q°q1q1)  + 6,*(PUP1P'VqV) 


On InO  n In^n 1 ' 


12  3 12  3 


12  3 12  3 


+ e.*(p1p1p°q°qcq1)  ♦ b r*(p0p°p1q1q1q0)  + 


12  3 12  3 


12  3 12  3 


6 *(p1p°p1q°q1q°) 

6 12  3 12  3 


+ 0 *(p°p1p1q1q°q°)  + 0 *(p1p1p1q°q°q°) 

7 123123  8 123123 

Since  *e4 >e6 »67*08  = 0 and  62’63,B5  = lf  the  above  equation 


1 


. A 
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reduces  to: 


P( ^ ) = plpOpOgOqlql  + pOplpOqlqOql  + pOpOplqiqlq 
V1/  1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 


= pqq  + p q q + p q q 

1 2 3 2 1 3 3 1 2 


Similarly, 


’(O'  pip< 


q + p p q + p p q 
3 1 3 2 2 3 1 


and 


’(0 ' p. 


p p 
2 3 


Thus,  the  probability  of  at  least  1 success  out  of  7>  events  is 
computed  by 


-(:)  • -(;)  • -0 


■pqq  + pqq  + pqq  + ppq  + ppq  + ppq  + ppp 

123  213  312  123  132  231  123 


In  instances  in  which  only  the  ait  least  case  is  required,  the 
solution  can  be  determined  with  the  formation  of  fewer  products,  thus 
reducing  the  number  of  computations.  This  procedure  can  be  utilized 
by  taking  advantage  of  complementary  logic. 


N 

It  follows  that  when  X < y , then  Eq.  (6)  should  be  used  to 

compute  the  solution.  In  the  preceding  case  for  N = 3 and  X = 1,  if 
the  only  requirement  had  been  the  at  least  case,  then  the  solution 
would  be  obtained  using  Eq.  (6).  Thus, 

i' ® 

i=  1 


1 - P 


(o) 


1 ' 

12  3 


The  abo-,e  procedure  requires  only  1 product  combination  rather 
than  the  previously  computed  7 combinations.  The  choice  of  using 
Eq.  (5)  or  Eq.  ^6)  can  best  be  illustrated  by  a more  realistic  ex- 
ample in  which  case  N = 18  and  X = 9.  If  Eq.  (5)  were  to  be  used. 


then 


18  / \ 
5 (")  ■ 


155,382  combinations  would  be  required  whereas  the 

8 


choice  of  Eq.  (6)  would  require  only 


£ (”)  • 
1=0 


= 106,762  combinations. 


A FORTRAN  computer  code,  NX,  has  been  written  for  efficiently 
computing  the  probability  of  at  least,  exactly  or  both  cases  as 
expressed  in  Eq.  (2).  This  code,  which  takes  advantage  of  the 
complementary  logic  condition  where  applicable,  is  given  in  Appendix 
A.  In  addition,  a similar  code  in  BASIC  is  also  given  in  Appendix  B 
for  implementation  on  mini -computers  utilizing  this  language. 


The  following  sample  problem  was  run  using  the  FORTRAN  code 


given  in  Appendix  A: 


Compute  the  probability  of  exactly  5,  exactly  6,  exactly  7 . . . 
exactly  13  successes  out  of  13  trials.  The  individual  probabilities 
p.  are  as  follows: 


p = 0.30,  p = 0.16,  p = 0.22,  p = 0.28, 

12  3 4 

p = 0.33,  p = 0.51,  p = 0.44,  p = 0.51, 

7 8 9 10 

p = 0.73,  p = 0.94  . 

12  13 

The  code  produced  the  following  solution: 


0.10,  p 

6 

= 0.68, 


Probability  of 

exactly 

5 

successes 

out 

of 

13 

trials 

= 

0.237 

Probability 

of 

exactly' 

6 

successes 

ouf 

of 

13 

trials 

= 

0.239 

Probability 

of 

exactly 

7 

successes 

cut 

of 

13 

trials 

= 

0.165 

Probability 

of 

exactly 

8 

successes 

out 

of 

13 

trials 

= 

0.079 

Probability  of  exactly 

9 

successes 

out 

of 

13 

trials 

- 

0.026 

Probability 

of 

exactly 

10 

successes 

out 

of 

13 

trials 

= 

0.006 

Probability 

of 

exactly 

11 

successes 

out 

of 

13 

trials 

= 

0.001 

Probability 

of 

exactly 

12 

successes 

out 

of 

13 

trials 

= 

0 . 000 

Probability 

of 

exactly 

13 

successes 

out 

of 

13 

trials 

= 

0 . 000 

Probability 

of 

at  least 

5 

successes 

out 

of 

13 

trials 

= 

0.754 

0.39, 


This  example  required  7,099  product  combinations  to  be  formed  and 
computed.  Run  times  were  0.6  second  and  2.0  seconds  on  a CDC  Cyber  76 
and  a CDC  Cyber  173,  respectively . Using  the  code  listing  in  Appendix 
B,  the  identical  results  were  obtained  from  a WANG  22C0-T  mini- 
computer with  a run  time  of  approximately  83  minutes  and  from  a WANG 
2200 -VP  with  a run  time  of  4.5  minutes. 

For  further  information  on  the  subject,  the  reader  is  referred  to 
Feller,  William , "An  Introduction  to  Probability  Theory  and  It '8 
Applications" , Nev  York,  NY:  John  Wiley  ■$  Sons,  1968. 


III.  CONCLUSIONS 

A general  expression  is  presented  for  computing  the  probability 
of  X kills  (successes)  out  of  N trials,  with  unequal  or  equal  indi- 
vidual probability  of  kill  (success)  values.  Computer  codes  in  both 
FORTRAN  and  BASIC  are  presented  which  yield  the  desired  solutions 
using  appropriate  main-frame  or  mini -computers . 
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Input  instructions  for  FORTRAN  version  of  NX  code 


Card  # 

Format  Type 

Description 

1 

215 

Number  of  trials  N,  number 
of  successes  X 

2 

16F5.1 

N probabilities  of  success 

3 

15 

Code  to  select  method  of 
computation; 

= 1 to  compute  exactly  case 

= 2 to  compute  both  at  least 
and  exactly  case 

= 3 to  compute  at  least  case 

14 


PROGRAM  NX  ( INPUT .OUTPUT ,TAPE5= INPUT ,TAPE6=OUTPUT ) 
DIMENSION  13(201),  P(200),  A(3),  Z(2) 

INTEGER  X,  R3 

1 Z ( 1 ) =5HEVENT 
IZ  = 1 

Z (2 ) =6HEVENTS 
A ( 1 ) =7HEXACTLY 
A (3 ) =8HAT  LEAST 
A (2 ) =10HBOTH  TYPES 
READ  (5,16)  M.X 
IF  (X.EQ.l)  GO  TO  2 
IZ=2 

2 IF  (M.GE.X)  GO  TO  3 
STOP 

3 READ  (5,17)  (P(I),I=1,M) 

N = 0 

DO  4 1 = 1, M 

IF  (P(I).EQ.O.)  GO  TO  4 
N=N  + 1 
P <N ) =P  ( I ) 

4 CONTINUE 

IF  (X.GT.N)  GO  TO  15 
T 1=M-X 

IF  (T1.NE.0.)  GO  TO  5 
Tl  = l . 

5 CONTINUE 

READ  (5,16)  R3 
GO  TO  (6,7,11),  R3 

6 N1=X 

CALL  D (N,N1,P,S> 

WRITE  (6,13)  X,Z(IZ),M,S 
STOP 

7 S2=0 . 

DO  9 N1=X,N 
CALL  D (N,N1,P,S) 

IF  (Nl.EQ.l)  GO  TO  8 
IZ=2 

8 WRITE  (6,19)  N1,Z(IZ),M,S 

9 S2=S2+S 

IF  (X.EQ.l)  GO  TO  10 
IZ=2 

10  WRITE  (6,20)  X , Z ( IZ ) ,M ,S2 
STOP 

11  S2=0 . 


15 


■ I 


IF  (2*X.LT.N)  GO  TO  13 
DO  IS  N1=X ,N 
CALL  D (N,N1,P,S) 

S2=S2+S 

12  WRITE  (6,20)  X,Z(IZ),M,S2 
STOP 

13  DO  14  NN1=1,X*N1=NN1-1 
CALL  D (N,N1,P,S) 

S2=S2+S 

14  CONTINUE 
SS2=1.-S2 

WRITE  (6,20)  X,Z(IZ) ,M,SS2 
STOP 

15  WRITE  (6,21)  X 
STOP 

16  FORMAT  (215) 

17  FORMAT  (16F5.1) 

18  FORMAT  (5X,22HPR0BABILITY  OF  EXACTLY , 15 ,A7 , 6HOUT  OF , 14 , 1H= ,F5 .3 ) 

19  FORMAT  ( 5X , 22HPRDBAB IL IT  Y OF  EXACTLY , 13 ,A7 , 6HOUT  OF , 15 , 1H= ,F5 .3 ) 

20  FORMAT  (5X  .23HPR0BABILITY  OF  AT  LEAST , 13 ,A7 ,6HOUT  OF  , 14 , 1H= ,F5 .3 > 

21  FORMAT  (55H  PROBABILITY  OF  THIS  OCCURING  IS  ZERO  BECAUSE  YOU  DO  N 
1 , 16HOT  HAVE  AT  LEAST ,F5 . 0 ,28H  NON-ZERO  PROBABILITY  VALUES) 

END 

SUBROUTINE  D (N,N1,P,S> 

DIMENSION  P (200 ) , 13(201) 

1 5 = 0. 

DO  2 1=1, N1 
13  < I ) =1 

2 CONT1NUL 
NN1=N1+1 

DO  3 I=NN1,N 
13 ( I > =0 

3 CONTINUE 

4 Al=l. 

DO  6 1=1, N 

IF  (13(1) .EQ.0)  GO  TO  5 
A1=A1*P ( I ) ' 

GO  TO  6 

5 A1=A1* ( 1 . -P ( 1 ) ) 

6 CONTINUE 
S=S+A1 

IF  (N.NE.N1)  GO  TO  7 
RETURN 

7 1 = 0 


a i=i+i 

IF  <I3(N+1-I) .GT.O)  GO  TO  9 
IF  (I.NE.N)  GO  TO  8 
RETURN 

9  IF  ll.LT.2)  GO  TO  10 
I4=N+1-I 
13 ( 14  ) =0 
14=14+1 
13(14) =1 
GO  TO  4 

10  I9=N 

IF  (13(19-1). LT.l)  GO  TO  13 
J = 1 

11  J=J+1 

IF  (I3(I9-J> .EQ.O)  GOTO  IS 
IF  (J.NE.N-1)  GO  TO  11 
STOP 

12  I9=I9-J 

IF  (J-N1.NE.0)  GO  TO  13 
RETURN 

13  17=19-1 
J = 0 

14  J=J+1 

IF  (I3(I9-J) .GT.O  i GO  TO  15 
IF  (J.NE.I7)  GO  TO  14 
RETURN 

15  I3=I9-J 
13 ( IS ) =0 
18=18+1 
13 ( 18 ) =1 

IF  (J.EG.l)  GO  TO  4 
J = 0 
N3  = 0 

16  J=J+1 

IF  (I3(18+J) .EQ.O)  GO  TO  17 
N3=N3+i 

17  IF  (I8+J-N.LT.0),  GO  TO  16 
L = 0 

18  L=L+1 

IF  ( 13 ( I8+L ) .NE . 0 > GO  TO  4 
13(18 +L)=1 
13 (N+l-L ) =0 
IF  (L.EG.N3)  GO  TO  4 


10  REM  NX 

SO  COM  I3(201),P(200),Ai(3)60,Z*(4) 

302  PROBABILITY  OF  EVENT  ##  IS 

401  PROBABILITY  OF  EXACTLY  ####  ######  OUT  OF  ####  = #.«##» 

502  PROBABILITY  OF  AT  LEAST  ####  #«##*«  OUT  OF  ####  = #.#### 

6 OX#########  ################  BEING  CALCULATED,  PLEASE  BE  PATIENT 
70  Z< < 1 ) = *EVENT “ : Z*t2)="EVENTS" : Z*(3)  = "COMBINATION  IS":  Z*(4)=*C0MB 
INAT  IONS  ARE".*  2 = 1 : STR < Z* <3 ) , 16 , 1 ) =HEX < 08 > 

80  PRINT  HEX (03) 

90  A* ( 1 ) = “EXACTLY " : A5(3)="AT  LEAST":  A*(2)="B0TH  TYPES" 

100  PRINT  "PROGRAM  IS  DESIGNED  FOR  A MAXIMUM  OF  200  EVENTS"  , HEX (0A0A ) 
110  INPUT  "TOTAL  NUMBER  OF  EVENTS", M:  IF  M<=0  THEN  110".  INPUT  "NUMB 
ER  OF  EVENTS  FOR  WHICH  YOU  WISH  TO  CALCULATE" ,X 
ISO  IF  X=1  THEN  130:  Z=2 

130  IF  M>=X  THEN  140:  PRINT  HEX  < OA ),  "ERROR  GOTO  110 

140  PRINT  HEX (04):  FOR  1=  1 TO  M 

150  PRINTUSING  30,1;  ."INPUT  P ( I > : IF  P(I)*<1-P(I  > >>=0  THEN  160:PRINT  HE 

X ( OC ),* ILLEGAL  VALUE,  TRY  AGAIN  ":  GOTO  150 

160  NEXT  I 

170  PRINT  HEX (03) 

180  N=0:  FOR  I = 1 TO  M:  IF  P(I)=0  THEN  190:  N=N+1 :P (N ) =P ( I ) 

190  NEXT  I:  IF  X>N  THEN  660 
200  DEFFN  0:  PRINT  HEX (03) 

E10  PRINT  "IF  YOU  WISH  TO  COMPUTE  THE  PROBABILITY  OF  EXACTLY  X EVENT 
S OUTOF  N EVENTS,  KEY  IN  A 1":  PRINT  HEX(OA) 

S20  PRINT  "If  YOU  WISH  TO  COMPUTE  THE  PROBABILITY  OF  AT  LEAST  X EVENT 
S OUTOF  ti  EVENTS  WITH  EXACT  INDIVIDUAL  PROBABILITIES,  KEY  IN  A S":  PRI 
NT  HEX  < OA ) 

230  PRINT  "IF  YOU  WISH  TO  COMPUTE  THE  PROBAB1LTY  OF  AT  LEAST  X EVENT 
S OUTOF  N EVENTS,  KEY  IN  A 3" 

240  FRINT  HEX  < OAOA ) 

250  INPUT  “YOUR  OPTION  NUMBER", R3:  PRINT  HEX (03) 

260  ON  R3  GOTO  270 ,300 ,400 :PRINT  "ILLEGAL  VALUE":  GOTO  210 

270  X1=X  : N 1 , N2  =N : GOSUB  '79:  PRINTUSING  60 ,C4 , Z* < Z1 ) : PRINT  : A=C4 

280  N1=X:  GOSUB  '99:  PRINTUSING  2°0 ,X , Z$ ( Z ) ,M ,S :STOP 

2902  PROBABILITY  OF  EXACTLY  ####  ######  OUT  OF  ####  = #.#### 

30c  a=o:  xi=x:  ni,n2=n:  gosub  '63:  printusing  6o,a,z*(Z1):  print 

310  S2  = 0:  FOR  Ni-x  TO  N:  GOSIJB  '99 

320  IF  Nl=l  THEN  330:  Z=2 

330  PRINTUSING  40  ,N1 , Z$ ( Z ) ,M ,S 

340  S2=S2*S:  NEXT  Ni:  S=0 

350  IF  N1=M  THEN  370:  N1=N1+1 

360  PRINTUSING  40  ,N1 , Z* ( Z ) ,M ,S : GOTO  350 

370  Z = 1 

380  PRINT  HEX(OA):  IF  X=1  THEN  390:  Z=2 
390  PRINTUSING  50  ,X , Z$ ( Z ) ,M ,S2 :STOP 

400  S2=0:  IF  2*X<=N  THEN  430:  X1=X:  N1,N2=N:  GOSUB  '69 
410  PRINTUSING  60 ,A , Z$ ( Zi ) :PRINT 
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420  FOR  N1=X  TO  N:  GOSUB  '99:  S2=S2+S:  NEXT  Nl:  PRINTUSING  50,X,Z*(Z> 
,M,S2:  STOP 

430  A,X1=0:  N2=X-l:Nl=N:  GOSUB  '69 

440  PRINTUSING  60  ,A , Z$ ( Z1 ) :PRINT 

450  IF  XOO  THEN  460:S2  = 0 :GOTO  470 

460  FOR  N1=0  TO  X-l:  GOSUB  '99:  S2=S2+S:  NEXT  Nl 

470  PRINTUSING  50 ,X , Z$ ( Z ) ,M , 1-S2 : STOP 

48 0 DEFFN'99:  S=0 

490  FOR  I = 1 TO  Nl:  I3(I)=l:  NEXT  I 
500  FOR  I = Nl  + 1 TO  N:  13 < I > =0 : NEXT  I 

510  Al=l:  FOR  I = 1 TO  N:  IF  13 ( I ) =0  THEN  520:  A1=A1*P<1):  GOTO  530 
520  A1=A1* ( 1-P ( I ) ) 

530  NEXT  I:  5=S+A1:  IF  NON1  THEN  540:  RETURN 
540  1=0 

550  1 = 1 + 1:  IF  13 (N  + l-I ) >0  THEN  560*.  IF  ION  THEN  550:  RETURN 

560  IF  I <2  THEN  570:  I4=N  + 1-I*.  I3(I4)=0:  14=14+1;  13(14)=l:  30TC  510 

570  19=N:  IF  13(19-1X1  THEN  600:  J = 1 

580  J=J  + l:  IF  13(19-1)00  THEN  580 

590  I9  = I9-J IF  J-N1O0  THEN  600:  RETURN 

600  17=19-1:  J=0 

610  J=J  + i:  IF  13 ( I9-J ) >0  THEN  620:  IF  JOI7  THEN  610:  RETURN 
620  I8  = I9-J : I3(I8>=0:  18  = 18  + 1*.  I3(I8)=i:  IF  J = 1 THEN  510:  J ,N3  = 0 
630  J=J  + 1*.  IF  13 ( 18 +J  ) =0  THEN  640:  N3=N3+1 
640  IF  I8+J-N<0  THEN  630:  L=0 

650  L=L  + 1!  IF  13  ( 18+L  K>0  THEN  510:  I3(18+L)=l:  I3(N  + 1-L)=0*.  IF  L=N3  T 
HEN  510:  GOTO  650 

66C  PRINT  ‘PROBABILITY  OF  THIS  OCCURING  IS  ZERO  BECAUSE  YOU  DO  NOT  HA 
VE  AT  LEAST*  *,X;*  NON-ZERO  PROBABILITY  VALUES  !**.  STOP  . 

670  DE^FN'69:  A = 0:  IF  XiOO  THEN  680*.  Xl,A  = l:  IF  X1<N2  THEN  680:  IF  N2 
<-*0  THEN  710:  A=A  +N : GOTO  710 

680  IF  XlOl  THEN  690:  A=A+N : Xl=2:  IF  X1<N1  THEN  690:  GOTO  710 
690  FOR  J=X1  TO  N2 : GOSUB  '79*.  A=A+C4:  X1=X1  + 1 
700  NEXT  J 

710  Zl=4:  IF  AO  1 THEN  720:  Zl=3 
720  RETURN 

73C  DEFFN'79*.  C1,C2,C3  = 1*.  IF  Nl<2  THEN  740:  FOR  1=2  TO  Nl:  C1=C1*I:  NE 

XT  I 

740  IF  N1-XK2  THEN  760 

750  FOR  1=2  TO  Nl-Xi:  C2=C2*I:  NEXT  I 

760  IF  Xl<2  THEN  780 

770  FOR  1=2  TO  XI:  C3=C3*I:  NEXT  I 

780  C4=C1/<C2*C3> 

790  Zl=4:  IF  C4<>  1 THEN  800:  Zl=3 
800  RETURN 
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U.  S.  ARMY  ARMA 
U.  S AR> 
ABERDEEN 


^RTMENT  OF  THE  ARMY 

RESEARCH  AND  DEVELOPMENT  COMMAND 

(.LISTIC  RESEARCH  LABORATORY 
ING  GROUND.  MARYLAND  21005 


DRDAR- BLV 


12  June  1980 


SUBJECT:  Memorandum  Report  ARBRL-MR-02890,  Jan  1979,  by  M.  B.  Danish 
and  J.  T.  Hundley 


SEE  DISTRIBUTION 


1.  The  subject  report  contains  two  computer  programs.  The  authors  have 
learned  of  an  algorithm  which  accomplishes  the  same  results,  but  much 
faster,  as  those  obtained  from  the  codes  given  in  subject  report.  The 
algorithm  was  published  in  a formal  report  at  the  Naval  Surface  Weapons 
Center,  Dahlgren,  VA  in  1975  but  distribution  was  at  that  facility 
only,  hence  the  information  of  importance  was  not  known  to  these  authors 
at  BRL  until  after  ARBRL-MR-02890  was  published.  It  should  be  noted 
that  the  original  algorithms  used  in  the  BRL  version  were  reported  on 

at  BRL  in  May  1974. 

2.  In  the  interest  of  computing  efficiency,  we  have  utilized  the  algorithm 
developed  by  M.  A.  Thomas  and  A.  E.  Taub  at  the  Naval  Surface  Weapons 
Center  and  prepared  a FORTRAN  and  BASIC  version  to  replace  those  versions 
originally  appearing  in  ARBRL-MR- 02890.  The  Thomas/Taub  algorithm  is 
extremely  efficient  and  can  be  considered  to  be  orders  of  magnitude 
faster  than  the  BRL  versions. 

3.  The  reference  report  is  NSWC/DL  TN-DK-25/75,  Binomial  Trials  with 
Variable  Probabilities,  August  1975,  by  M.  A.  Thomas  and  A.  E.  Taub. 
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BRADSHAW  F.  ARMEN DT 

Chief,  Ground  Mobility  8 Firepower  Branch 
Vulnerability/Lethality  Division 


PROGRAM  NX ( INPUT* OUTPUT *TAPE5= I NPUT.TAPF6*0UTP0T) 
COMMON  T ( 20  *21 ) * H (20) 

PE AO  (5*900) K1  »N9 

IF (Kl.LE.O)GOTC  330 

WEAO  (5*901) (R(I9) ,I9sl*Kl) 

T(l*l)*l-«(1) 

T ( 1 *2) *R  ( 1 ) 

00  230  1 9*2 • K 1 
K9»I9>1 

00  220  J9*l*K9 
IF  (J9.NE.1)  GOTO  180 
T(I9.1)»T(I9-1,1)*(1-R(I9) ) 

GOTO  220 

1 80  IF(I9.EU.K9)  GOTO  210 

T(I9.J9)»T(I9-1»J9-1)*9(I9)*T( 19-1  * J9) * ( 1-R ( 19) ) 
GOTO  220 

210  T(I9,J9)aT(I9-l*J9-l)*R(J9) 

220  CONTINUE 

230  CONTINUE 
I9sK1 
55*0 
I C 1 *N9* 1 
IC2*K 1*1 

DO  300  J9=IC1»IC2 
WRITE (6 *902) J9-1.T(I9*J9) 

S5*S5*T ( 1 9. J9) 

300  CONTINUE 

WRITE (6*903) N9«S5 

330  STOP 

900  F ORMAT (2110) 

901  FORMAT(10F10.5) 

902  FORMAT ( • • * • PROB  OF  EXACTLY  * • I<**5X.F1 1 .7) 

903  FORMAT(*0»**PPO«  OF  AT  LEAST  ' * I* ♦ 1 OX *F 1 0 .6 ) 

ENO 


u 

1 m 


■Han 


10  COM  R(E0>:  DIM  T<30,31) 

SO  PRINT  HEX (03) 

30  INPUT  "TOTAL  NUMBER  OF  EVENTS" ,K1:  IF  K1<=0  THEN  SO 
40  INPUT  "NUMBER  OF  EVENTS  FOR  WHICH  YOU  WISH  TO  CALCULATE " ,N9 
50  PRINT  HEX ( OA ) 

60  FOR  19*1  TO  Kl: PRINT  "P ( " ; 19; " ) “ ; : INPUT  R< 19): NEXT  19 
70  PRINT  HEX (03) 

80  PRINT  "PLEASE  BE  PATIENT  AS  I AM  COMPUTING  FOR  YOU" 

90  T ( 1 , 1 ) *1  -R  ( 1 ) 

100  T ( 1 ,S ) *R ( 1 ) 

110  FOR  I9*S  TO  Kl 
ISO  K9=I9+1 
130  FOR  J9»l  TO  K9 
140  IF  J9<>1  THEN  160 

150  T ( 19 , 1 ) *T ( 19-1 ,1 ) * < 1 -R ( 19 ) ) : GOTO  E00 
160  IF  I9=K9  THEN  190 

170  T < 19 , J9 ) *T (19-1 , J9-1 > *R ( 19 )+T < 19-1 , J9>* ( 1-R< 19) ) 

180  GOTO  300 

190  T < 19 , J9 ) =T ( 19-1 , J9-1 ) *R  < J9 ) 

300  NEXT  J9 
S10  NEXT  19 
EEO  I9=K1 

S30  PRINT  HEX (03):  S5*0 

S40  FOR  J9*  N9+1  TO  Kl+1 

350  PRINTUSING  370,  J9-1 ,T ( 19 , J9 ) 

360  S5=S5+T ( 19 , J9 ) 

370  % PROB  OF  EXACTLY  ####  * -#*#.##*####  -###.###### 

380  NEXT  J9 
390  PRINT  HEX  < OA ) 

300  PRINTUSING  330,N9,S5 
310  STOP 

330%  PROB  OF  AT  LEAST  ##«  ■ #.####### 


